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SUMMARY 


This report describes a technique which combines theory and 
experiments for determining relaxation times in gases. The technique is 
based on the measurement of shapes of the bow shock waves of low-fineness- 
ratio cones fired from high-velocity guns. The theory presented in the 
report provides a means by which shadowgraph data showing the bow waves 
can be analyzed so as to furnish effective relaxation times. 


Relaxation times in air were obtained by this technique and the 
results have been compared with values estimated from shock tube measure- 
ments in pure oxygen and nitrogen. The tests were made at velocities 
ranging from 4600 to 12,000 feet per second, corresponding to equilibrium 
temperatures from 3590° R (1990° K) to 6200° R (3440° K) , under which 
conditions, at all but the highest temperatures, the effective relaxation 
times were determined primarily by the relaxation time for oxygen and 
nitrogen vibrations . 


INTRODUCTION 


The effects of molecular vibration and dissociation upon the 
thermodynamic properties of air (and its constituent gases) at elevated 
temperatures have been extensively studied, and accurate tabulations of 
the properties of air have been made at temperatures up to 15,000° K, 
assuming that the air is in chemical and thermodynamic equilibrium 
(see refs. 1, 2, and 3)* However, when the gas density is low and local 
velocities are high, or when temperatures vary rapidly along streamlines, 
the assumption of equilibrium can not necessarily be made, and the rate 
of adjustment of the gas properties must be taken into account. 

It is known that when air is heated suddenly in strong shock waves, 
the translational and rotational degrees of freedom of the molecules 
adjust to the new state with essentially no time lag, and that the 
molecular vibration and dissociation adjust considerably more slowly 
(ref. 4). When the latter type of adjustments require a time that is of 
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the same order as the time required for the gas particles to traverse a 
flow field being studied, it may be important to know accurately the ^ 

magnitude of this time, since the flow can be significantly influenced 
by it. 

Theories for calculating the variations in gas properties in 
nonequilibrium flow through normal shock waves have been investigated by 
various authors: the relaxation of molecular vibrations in oxygen and 

nitrogen in reference 5; and oxygen dissociation relaxation in references 
6 , and 8 . Reference 5 indicates that the mechanism of the relaxation 
process for vibrational excitation in a pure gas can be described theo- 
retically. Even for this comparatively uncomplicated process, however, 
the existing theory does not provide complete information for calculating 
numerical values of relaxation time, and experimental results are needed 
to evaluate collision cross-section constants appearing in the theory. 

At higher enthalpy, when the relaxation involves chemical reactions, and 
when the gas is a mixture, as in the case of air, the theory is consider- 
ably more complex, and is not generally sufficiently precise for use in 
solving flow problems, where accuracy is essential. Consequently, it is 
desirable to have available experimental results indicating relaxation 
times for a wide range of enthalpies and gas densities, for use both in 
solving engineering problems in gas flow and in evaluating the theory. 

The technique that has been employed in obtaining practically all 
existing relaxation time data at high enthalpy has been to measure in 
shock tubes the variations with distance or time of the local gas prop- 
erties in the region just behind ^the advancing shock wave (see refs. 

9, and 10) . 

This report describes another technique for obtaining relaxation 
times from experiments which, for certain test conditions, may be more 
convenient than the shock tube tests. The technique consists in firing 
from high velocity guns cone-shaped models of low-fineness-ratio and 
recording the shape of the bow shock wave. Effective relaxation times 
are then calculated from the wave shape data using a theory developed by 
D. R. Chapman that is presented in this report. 

Tests were conducted to obtain relaxation times in air. At the 
model velocities of the tests, the enthalpy behind the bow shock waves 
was such that the nonequilibrium effects were due primarily to the 
excitation of oxygen and nitrogen vibrations. Relaxation times for 
vibrations are available from shock-tube tests for these two gases 
(ref. 5). In the present report, the results from the shock-tube tests 
in pure oxygen and pure nitrogen were used to estimate effective relaxa- 
tion times in air, and these times were compared with those determined 
from the tests described herein. 

A limited amount of data was obtained in which small effects of 
oxygen dissociation could be observed on the bow wave shape. However, 
such data were not sufficiently precise to permit the calculation of 
relaxation times for the dissociation. 
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NOTATION 


C P 

pressure coefficient 

h 

enthalpy 

M 

Mach number 

P 

pressure 

P 

V 

radial coordinate parameter. — — 

tu g 

X 

> r 

radial coordinate 

• r b 

slant length of cone 

R 

Reynolds number, also gas constant for unit mass 

R b 

model radius 

t 

T 

temperature ratio across shock wave, — also time 

T i 

T 

absolute temperature 

T 

mean temperature, between frozen and equilibrium values 

U i 

free-stream velocity 

U 2 

resultant velocity in region between solid cone and shock wave 

u o 

average velocity behind shock wave 

u,v 

velocity components behind shock wave 

v m 

maximum possible velocity, expansion flow 

W 

resultant velocity just behind the shock wave 

y 

distance between actual shock wave and reference shock wave 
(see fig. 2(b)) 

y 

distance from solid cone surface (see sketch (a)) 

yb 

distance between shock waves for frozen and equilibrium flow at 
the distance r b along model cone (see fig. 2(h)) 

z 

compressibility factor 


h 


t 


a streamline angle behind shock wave 

h„/h 

3 enthalpy ratio divided by temperature ratio, — j — 

T 2/ T l 

7 ratio of specific heats for ideal gas 

7 e effective ratio of specific heats (see eq. (A2l)) 

S s yaw angle of solid cone 

yaw angle of bow shock wave 
S* boundary-layer displacement thickness 

0 V angular coordinate of bow shock wave 

0 e bow shock wave angle for equilibrium flow 

Qf bow shock wave angle for frozen flow 

0 S solid or model cone half angle 

AO difference between shock wave angle and solid cone angle > 0v - e s 

p density 

o increase in apparent shock wave angle due to projection, 0p - 0 W 

r relaxation time 

Subscripts 

e equilibrium 

E effective 

f frozen flow 

g oxygen in equilibrium, nitrogen frozen 

n normal shock wave 

N nitrogen 


0 


oxygen 
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p uncorrected for projection effects 

s model cone surface 

S constant entropy 

t parallel to shock wave 

w bow shock wave 

0 standard air conditions (l atm., 59° F) , also reference conical 

flow 

1 free stream 

2 behind the bow shock wave 


METHOD OF EXPERIMENT AND ANALYSIS 


The excitation of the internal degrees of freedom in a gas with high 
total enthalpy is evident as an increase in the specific heat and, when 
dissociation occurs, also as an increase in the compressibility. The 
effect of the increased specific heat is to decrease the temperature in 
the flow behind a shock wave compared with that of an inert gas at the 
same enthalpy. The density then increases approximately inversely with 
the temperature, while the pressure in most cases of flow through shock 
waves is relatively unaffected by the excitation of molecular vibrations 
and by dissociation. This increased density is opposed to a small extent 
by the effect of the increased compressibility accompanying dissociation. 

In the case of a two-dimensional oblique shock flow, in order to 
satisfy continuity requirements, the increased density of real air behind 
the shock wave results in a reduced total stream tube cross section 
between the shock wave and the solid surface, so that for the same flow 
turning angle, the shock wave must assume a more acute angle with the 
flow direction. The same is true for a cone. The geometry of the bow 
wave on a cone can thus be used as a measure of the state Qf the air 
behind the shock wave. 

A cone of large included angle is a particularly advantageous model 
for studying variations in gas properties, because the flow on a cone 
with an attached shock wave is subject to considerably simpler and more 
exact analysis than that on most other shapes adaptable to ballistic 
types of tests, such as blunt bodies with detached shock waves. 

A theory which relates the shock wave shape to the thermodynamic 
properties of the gas is given for cones of large included angle in the 
following section. It shows that in the case of such cones (selected so 
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that the flow is everywhere supersonic) the how shock wave will he 
conical if the gas is inert or is in equilibrium, and will he curved hy 
a predictable amount when nonequilibrium flow exists. 

To provide the experimental shock waves needed for this investigation, 
conical models of 52.5° and 55° half angle were launched from guns at 
speeds up to 12,000 feet per second through still air at pressures selected 
to give frozen, nonequilibrium, and equilibrium flows. The experimental 
technique is further described in a later section. 


Theory 


The theory on which the analysis of the tests is based was developed 
at the Ames Research Center by Dean R. Chapman. Since it has not been 
published elsewhere, it is presented in this report. This presentation 
can conveniently be divided into three parts: (l) the calculation of the 

bow shock wave cone angle for frozen flow, that is, for flow in -which the 
relaxation time is so long that none of the inert degrees of freedom 
(vibration and dissociation, in this case) are excited, so that the air 
behaves like a perfect diatomic gas; (2) the calculation of the bow wave 
angle for flow with the air in equilibrium behind the shock wave, corre- 
sponding to a relaxation time of zero; and (3) the calculation of the bow 
wave shape for nonequilibrium flow with various relaxation times between 
zero and infinity. 

Frozen flow .- Numerical solution of the Taylor-Maccoll equations for 
flow on an infinite cone furnishes an accurate means of calculating the 
geometry of the bow wave, as well as the various details of the flow field 
between the bow wave and the solid cone for frozen flow. In reference 11, 
such solutions are tabulated for cone angles up to 50°. The approximate 
method described in appendix A was used to obtain solutions for the 
slightly larger cone angles required in the studies which are the subject 
of this report. Figure 1 shows a comparison of the bow wave angles as 
functions of Mach number calculated for frozen flow by this approximate 
method and by the exact solution for two infinite solid cones, one having 
a half cone angle of 50° and on of 55° • In this figure the ordinate A0 
is the difference between the shock wave semicone angle, 0 W , and the solid 
surface semicone angle, 0 S . It is evident that the bow wave angles 
obtained in the approximate solution are very nearly the same as those 
found from the exact calculations . 

The exact solution for the 50° cone was obtained from reference 11 
in which 7 is assumed to be 1.405., whereas for the 55° cone the exact 
solution was carried out (by Chapman) for a 7 of 1.400. The same values 
of 7 were used in the approximate solutions to make these solutions 
comparable with the exact solutions. 
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Equilibrium flow .- It is apparent from reference 11 that the 
numerical solution for the flow on cones is lengthy and time-consuming 
even when perfect gas relations are employed. If the actual thermodynamic 
properties of air were included, solving the equations numerically would 
be considerably more difficult 1 and the solutions would lack generality, 
so that at any given Mach number, separate solutions would be required 
for each combination of free-stream temperature and pressure. As a more 
practicable approach to the problem of calculating the bow wave cone 
angles and some of the details of the flow behind the bow wave. Chapman 
investigated the approximate solution presented below. (Two other 
solutions with similar objectives have been described in refs. 13 and l4.) 
The solution is presented in terms of an "effective" ratio of specific 
heats (7-g) which is related explicitly to the thermodynamic properties of 

the gas. The solid cones to which the calculations apply are chosen to 
be as blunt as possible, consistent with the requirement that the bow 
wave be attached and have a conical shape in equilibrium flow. This 
bluntness leads to a flow which can be described to a high degree of 
accuracy considerably more simply than the general cone flow. The 
assumption is made that the temperature, density, and pressure behind the 
shock wave are constant in the entire region ahead of the solid cone. 

Expressions derived in appendix A relate the density ratio, the bow 
shock wave angle, and the pressure coefficient behind the shock wave for 
the flow of a hypothetical gas having a ratio of specific heats that is 
constant across the shock wave and equal to y^. The expressions that 

relate these variables are the same as those for a perfect gas. Because 
of the way that 7 is defined, at any given temperature ratio t, the 

.Cj 

density ratio p^p determined by equation (A22) for the hypothetical 
flow is the same as that for the actual flow of the general gas given by 
equation (A20) . Equation (Al6), giving the shock wave angle, and equa- 
tion (A23), giving the pressure coefficient, are the same for the actual 
flow and the hypothetical flow, and so, since the density ratios are the 
same, 0 W and Cp are the same in the two flows (for the same solid cone 
angles and temperature ratios) . However, the enthalpy ratios are not the 
same, and as a result, the Mach numbers of the two flows are different 
and related by equation (A26) . 

As the first step in solving specific conical flow problems for the 
general gas, the solution is obtained for the mathematically simpler 
hypothetical flow, in terms of the effective y and the effective Mach 
number. The solutions corresponding to semicone angles (0 S ) of 52*5° 
and 55° are shown in figures 3 through ]_. With these solutions in 

■''Subsequent to the derivation of relations and the evaluation of the 
parameters for the cone flows pertaining to this report, reference 12 
became available describing the numerical integration of the Taylor- 
Maccoll equation for equilibrium flow which was programmed on an IBM 70^- 
computer. Such a solution, incorporating the proper free-stream conditions 
and solid cone angles, had it been available could have been used to 
obtain the cone-flow parameters required for the studies discussed below. 
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graphical form, as illustrated, for any assumed conditions which define 
the states of the gas on each side of the shock wave, 7^ is calculated 

from equation (A2l), and values of M , 0 W - 0 S , and Cp are read from 

1 E 

the graphs. Then, from equation (A26) the free-stream Mach number of the 
actual flow is obtained, and all of the required quantities have been 
evaluated for the equilibrium flow. 

Appendix B illustrates the calculation of the bow shock wave angle 
and pressure ratio for a 52.5° solid cone at a free-stream temperature of 
5^0° R at one Mach number and pressure. Results of this and similar 
calculations for other Mach numbers and pressures are plotted in figure 8 
which shows A0 , the difference between the bow shock wave semicone angle 
and the solid body semicone angle, as a function of Mach number. When 
oxygen dissociation occurs, the equilibrium bow wave angle is a function 
of the pressure, as indicated by the lower curves. The latter curves axe 
for constant pressure behind the shock wave. 


Nonequilibrium flow .- In this section an expression is derived from 
which the coordinates of the bow shock wave on a cone can be calculated 
for nonequilibrium flow. This flow has the characteristic that the time 
required for the gas to attain a density near equilibrium, after passing 
through the shock wave, is of the same order as that required for the gas 
particles to move a distance typical of the dimensions of the flow field 
being studied. 

In the following development, which first considers an infinite cone, 
three assumptions are made: (a) the flow properties (p, u, h) between the 

bow wave and cone surface differ by 
only a small amount from correspond- 
ing properties in some reference 
conical flow (p Q , u Q , h Q ); (b) this 
reference conical flow (which could 
be, e.g., either the frozen conical 
flow or the equilibrium conical flow) 
corresponds to a cone angle 0 S 
sufficiently blunt that all proper- 
ties p G , u Q , h Q are essentially 
constant between the bow wave and 
cone surface; and (c) the free-stream 
Mach number is sufficiently high that 
0 W - 0 S is small compared to 0 S 
(so that sin 0 = sin 0 W ) . Referring 
to the stream tube in sketch (a), 
from the equality of the free-stream 
mass flow in, 2jtp 1 u 1 r w sin 2 0 w dr w , 
and the flow out, 2rtpur sin 0 dy, 
and noting that pu may be expressed as a function of r w , 
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r sin 0 dy = 
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( 1 ) 
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i u i 

pu 


By applying this equation first to the nonequilibrium flow, then to the 
reference conical flow, and subtracting (noting from assumption (c) that 
sin 0 = sin 0 W = sin 0 Wq ) we have 


p o u o r 
p 1 u 1 sin 0 Wq 





Pq u o ~ pu 

pu 


T \ft r \r 


(2) 


For a particle flowing along a streamline the flow properties are regarded 
as a function of the enthalpy. From assumption (a) of small departures 
from a reference conical flow. 


pu 


Po u o + 


d(pu) 

dh 


(h 


h Q ) 


( 3 ) 
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equation (2) becomes, upon disregarding second-order terms in the inte- 
grand. 


PQ U Q 

Pl U l 


dh r 

^(pu)J o sin 0 w o 





h ~ h o 
Po u o 




The definition of relaxation time r is taken as 


w 


dh ^e ~ h 
dt t 


Since dh/dt = u dh/dr along a streamline, this equation for constant r 
and constant u = u Q (assumptions (b) and (a)) may be integrated between 

r = r TT and r = r to obtain 
w 

h e - h = (h e - h v )e"( r " r w)/^o (6) 


This equation shows that at r = oo, h = h e , as should be the case, and 
at r = r w , h = h w (by definition of hy.) . By noting that 
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h h Q - h e h Q (h 0 h) 


r - r. 


w 


= h e - h Q - (h e - h w )e Tu o 


It follows from equation (4) that, with x = and P = 


ru^. 


e - e 0 * ^ ~ y "° 

O r 


r 1 

-* I 




( 7 ) 


(8) 


Applying this now to the case of a nonequilibrium flow near a reference 
frozen flow, h Q = h-f, said h w will everywhere be nearly hf so that 
equation (8) integrates to 

e - e f = CF(P) ( 9 ) 


where C 


and 


is a constant equal to i 


(he ~ h f ) PjUjSin 0 f 


( Y 

' 8h 


_8(pu)_ 

r p + p - 1 ) 



( 10 ) 


In the limit of frozen flow T -+ °°, P -► 0, and F(P) -» (l/3)P, which gives 
in the limit, F(0) = 0, thus satisfying equation ( 9 ) • We have 


0 ~ 0 f _ CF(P) 
e e ~ e f d e ~ 9 f 


(ID 


p 

so that in the opposite limit of equilibrium flow, t -*■ 0, P = — * 00 , 

1u o 

and Q must approach 6 e . We use this requirement to evaluate C from 
the above equation (noting that F(oo) = l) . 


1 = 


C 


0 e - 0f 


( 12 ) 
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consequently , 


J- ^ = F(P) = 1 - ■% (e“ P + P - 1) 

9 e ~ e f P 


( 13 ) 


or 



(i*0 


k 

2 

2 
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This furnishes an expression from which the coordinates of the how shock 
wave on an infinite cone can he calculated when the gas behind the shock 
wave is not in equilibrium. 


" ®e 2 

0 f - 0 e = ( r/ tuq ) 2 


(' 


e 


-r/ tu 0 



) 


(15) 


In this equation 0 W is the angular coordinate and r is the radial 
distance to points on the shock wave. Also, u 0 is the average velocity- 
in the flow behind the shock wave, and t is the relaxation time. In 
the derivation presented, t is the time required for the enthalpy of the 
gas to reach a specified fraction of its total variation from frozen to 
equilibrium. If t is assumed to be a constant during the relaxation, 
this fraction corresponds to the state when the value of the state param- 
eter deviates by l/e times the initial deviation from the equilibrium 
value . 

It is convenient to alter the form of equation (15) slightly in 
applying it to the case of finite cones. In this altered form, the 
distance ratio y/y-^ is calculated as 


As indicated in figure 2(b), y is the distance measured from the equili- 
brium shock wave to the curved shock wave at radial distance r, and y^ 
is the distance between the frozen and equilibrium shock waves at radial 
distance r^. This distance ratio y/y^ has been calculated and is shown 

graphically in figure 9 as a function of the radial coordinate r/r^ for 
values of the reciprocal relaxation distance ratio r^/m from zero to 

infinity. At the extreme values of this ratio, zero and infinity, the 
solutions coincide, respectively, with the solutions for frozen and 
equilibrium flow behind the shock wave. 


^b 


y 


( 
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Experiments 


The theory presented in the preceding sections shows that, from 
information on the included angle of the how shock wave of low- fineness- 
ratio cones, the state of the gas behind the shock wave can be deduced, 
and in cases where the relaxation adjustments in the gas occur in a 
portion of the flow which is comparable in size with the region of flow 
on the conical portion of the model, the relaxation time in the gas can 
be estimated from the contour of the (curved) shock wave. Experiments 
were undertaken in which the state parameters of the gas would vary 
appreciably between the equilibrium and frozen states, and test conditions 
were selected for which it was expected that nonequilibrium flow would 
occur. The experiments consisted of launching conical models (of 52.5° 
and 55° half angle) in free flight at high velocity in the Supersonic 
Free-Flight Wind Tunnel. The following information was required from the 
experiments: (l) well-defined pictures of the bow wave from which the 

wave shape or wave angle could be measured and (2) measurements of the 
free-stream air properties and the total velocity with which the measured 
wave geometry could be correlated. 

Apparatus and models .- The pictures of the bow waves were provided 
by spark shadowgraphs which are obtained at each of nine instrumented 
stations in the wind tunnel test section (see ref. 15). For the tests 
reported, three smooth-bore guns were used to launch the models: a 
1-3/^- inch-bore powder gun, a 37 ~™i light-gas gun, and a caliber 50 light- 
gas gun. The latter two used helium (compressed in a shock- tube reservoir) 
as the propellant, as described in reference 15. 

The models, examples of which are shown in figure 10, were short 
blunt cone-cylinders, the forward portions of which were machined from 
aluminum alloy and were bonded to the cylindrical nylon rear portion with 
an epoxy resin. (it appeared that screw fastenings between the parts 
would often promote model breakage in the gun, as a result of stress 
concentrations.) The flared afterportions of two of the models in 
figure 10 (which were for use in the light-gas guns) sheared off before 
the models left the launch tubes. 

The models were accurately measured and inspected with a contour 
projector (at 20x magnification), and only those models were used which 
had exact straight-sided cones without bluntness or concavity. Slight 
variations of the included cone angle from the nominal angle were noted, 
amounting to a maximum of 0.2°, but because of the way in which the data 
were analyzed, these variations had no effect on the results of the tests. 

Tests.- The majority of the data presented in this report has been 
obtained by firing the models into still air at pressures ranging from 
1 atmosphere to 0.057 atmosphere, and at velocities between 7 ,260 
12,000 feet per second. In these tests the static temperature of the air 
was about 530° R. 
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Tests were also conducted by firing into a supersonic air stream of 
Mach number 3 in the Supersonic Free-Flight Wind Tunnel, which is 
described in reference 15 . In these tests, for which the air stream 
static pressure was approximately 0.10 atmosphere and the static tempera- 
ture was approximately 190° R, data were obtained at Mach numbers from 
6.79 to 13*55> corresponding to equilibrium temperatures behind the bow 
shock wave ranging from 1510° R to 4460° R. 

A relatively large number of rounds were fired in the test program 
but the amount of data obtained was not large, primarily because of two 
problems: breakage of the model in the launch tube, and oscillations of 

the model attitude in flight. The oscillations were particularly trouble- 
some when the models were fired into a partial vacuum. In this case the 
gun muzzle blast often produced a large disturbance to the model attitude, 
and periods of oscillation were so long that at the instants when the 
shadowgraphs were recorded, the probability that the model would be at a 
small angle of attack was low. 

Reduction of data .- The data obtained in the tests consisted of the 
time-distance histories of the model flights, together with spark shadow- 
graphs at the nine instrumented stations. With (the Mach number 3) air 
flow in the test section, the free-stream conditions were determined from 
measured temperatures and pressures in the settling chamber in conjunction 
with the wind-tunnel calibration. 

At each station the angles of yaw and pitch of the models were 
measured from the shadowgraphs and plotted to reproduce the pitch and yaw 
histories of the flights. These plots were inspected in order to select 
the shadowgraphs in which the model was at a sufficiently low angle so 
that the measured bow wave shapes would be very nearly the same as those 
for the model at exactly zero angle of attack. As a guide to the maximum 
permissible angles of attack, the effect of yaw on a conical bow wave was 
considered. According to references 16 and 17, the bow wave of an infinite 
cone in flow with constant 7 (7 = 1 . 405 ) remains a circular cone with 

wave angles unchanged even at relatively large angles of yaw. Only a 
second-order solution indicates a small eccentricity to the circular cone 
wave shape. Reference l6 also shows that when the solid cone is yawed, 
the axis of the bow wave cone is generally no longer coincident with the 
solid cone axis. This is illustrated in figure 11 where the ratios of 
the yaw angles of these axes are plotted as functions of the solid cone 
semi vert ex angle 0 S . 

When yawing of the model is moderate and occurs only in the plane 
parallel to one shadowgraph film plate, the included angle of the bow 
wave in the shadowgraph is (according to the solution in which 7 is 
constant) the same as that for zero yaw angle. 

In addition to the possibility of an aerodynamic effect of yaw, an 
optical effect is present such that in a plane normal to the plane of 
yawing, a shadowgraph plate will record the projection on that plane of 
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the bow wave which has an apparent angle 0-p larger than 9w according 
to the relation 


cos 5 "J cot 2 0 w - tan 2 S 

■ban 0 -^ = - 

p sin 2 S - cos 2 & cot 2 0„ 

■where 6 is the yaw angle of the how wave. The difference between the 
actual shock- wave cone angle, 0^ and the projected angle, 0p, a = 0p - 0 V , 

is plotted in figure 12. By the use of figure 11, extrapolated to include 
model cone angles of 52 . 5 ° 55° and figure 12, it was generally possible 

to select shadowgraphs from most of the test flights in which the effect 
of this type of distortion would be less than 0.10°. A correction to the 
data was applied when the effect of this distortion was significant. 

Spark shadowgraphs which showed well-defined bow waves and met the 
requirements that the model yaw was small were selected and accurately 
measured. Typical shadowgraphs are shown in figure 13. The coordinates 
of a series of points along the front of the bow wave image were read 
parallel and perpendicular to the bisector of the bow wave. A film- reading 
machine fitted with a microscope was used to measure the coordinates, 
which could be read to an accuracy of 0.001 inch in the streamwise direc- 
tion and to 0.005 inch laterally. Figure l4(a) shows an example of the 
bow wave contour that was read from a shadowgraph in this manner and 
plotted. 

The bow shock waves in the shadowgraphs were not always symmetrical 
about a bisecting axis, but instead had different curvatures on the two 
branches. Sometimes the reasons for this asymmetry were not evident, but 
in other cases it could be attributed qualitatively to the effect of a 
small yaw angle, which would, in nonequilibrium flow, produce different 
rates of excitation of the vibrations and dissociation on the two flow 
regions, the higher temperatures and pressures on the windward side causing 
a decrease in the relaxation time. Because the theory that has been devel- 
oped required the assumption that the wave angles are symmetrical, when 
the bow wave was not symmetrical the curve was oriented symmetrically 
about the outer portions of the curve (just inside the estimated position 
where the model shoulder would influence the flow) . 

To calculate the relaxation time from the bow-wave data, the 
procedure described below was employed. This calculation is slightly 
indirect, but it tends to smooth errors in the location of individual 
points on the bow wave. Coordinates of the points read from the original 
shadowgraphs were transferred by simple rotation of the axes, so that the 
outer portion of the bow wave would be symmetrical about the horizontal 
axis. Polar coordinates of each point were calculated based on a location 
of the origin at the apex of the bow wave cone. An exact determination 
of the location of this apex sometimes proved to be the principal factor 
in limiting the accuracy of the final results. Although the models were 
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all examined carefully before firing to see that they were truly conical 
and not blunted or excessively pointed, some of the shadowgraphs evidenced 
local distortion, refraction effects, fuzziness or other causes of indis- 
tinct definition near the apex which would introduce errors. When it 
appeared that such an error might be large, the shadowgraph involved was 
not used. 

The polar coordinates of the points on the bow wave then were plotted 
as Qy as a function of radial distance, r/r-^. This type of plot, shown 

in figure l4(b) , was generally erratic near the origin, at r/r-^ less 

than about 0.3, hut became relatively smooth and consistent as the radial 
coordinates increased to unity. A curve was faired through the portion 
of the plotted points that defined a consistent relation, treating the 
upper and lower branches of the bow wave separately. From these two 
curves, a single mean curve was drawn to represent the best symmetrical 
average bow wave contour. At two or three values of r/r^ between 0.6 

and 1.0, the polar angle was read from this single curve and y/y^ 

was computed from the relation 


This value of y/y^, together with the corresponding coordinate r/r^, 
locates a point in the graph, figure 9> which determines a value of the 
parameter r^/ tu q • The relaxation time then is 


In this relation u 0 is the average particle velocity along a streamline, 
starting at the shock wave and extending downstream to a point where the 
gas is close to equilibrium or is leaving the region being considered. 

The calculation of an accurate value of u 0 is complicated by the 
necessity of knowing the details of local variations of the state of the 
air which in turn are functions of relaxation time. Preliminary calcula- 
tions indicated that the value of u Q obtained assuming that the air 
reached equilibrium with no time lag would not be greatly different from 
the value obtained assuming frozen flow with a bow wave angle equal to 
that for the equilibrium flow. Thus the calculation for the case of 
frozen flow would correspond to a solution for a more slender solid cone 
than that of the actual model, but for the same free-stream Mach number. 

It was estimated that the value of u Q determined in this manner would 
be accurate within 5 percent. From this it appeared that the degree of 
approximation involved in using this value of u Q would be well within 
the accuracy of the measurements for determining values of r; therefore, 
this value of u Q was used in reducing all of the data. 
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Boundary- layer effects .- The possible effects on the bow shock wave 
shape of the boundary layer on the conical model surface were investigated. 
In appendix C results of calculations are presented and discussed which 
show the maximum magnitude of the boundary-layer displacement thickness 
and the estimated effect on the bow wave angle in conical flow. On the 
basis of these estimates, it was concluded that the boundary layer had 
only a very small effect on the bow wave shape in practically all of the 
tests from which data were derived, and so no correction for this effect 
was made. 


RESULTS AND DISCUSSION 


Figures 15(a) through 15(h) show the experimental data that were 
used in the calculation of the relaxation times. In these figures, as 
in figure 8, the ordinate A0 is the difference between the angular 
(polar) coordinate of points on the shock wave and the half angle of the 
solid cone. Each pair of symbols joined by a solid line corresponds to 
a separate shadowgraph. When the two symbols are not coincident, the 
bow wave had curvature and the wave angle was measured on the shock wave 
at locations ranging from a point near the apex to a point near the 
maximum model radius. A total of 20 model shots is represented by the 
data shown in figure 15 . 

In addition to the experimental data, figure 15 shows (as in fig. 8) 
the curves calculated from the theory for frozen flow, flow with the 
oxygen vibrations fully excited, flow with dissociation frozen, oxygen 
and nitrogen vibrations fully excited, and flow with vibrations and 
dissociation in equilibrium. 

Figure 15 shows that most of the data lies in the region between the 
curves for frozen and for fully excited equilibrium flow. From this it 
is inferred that a major portion of the data contains regions of non- 
equilibrium flow. Since the data shown were obtained at various free- 
stream pressures, and since the relaxation time is a function of the 
local pressure (and the local pressure is a function of the stream pres- 
sure), it would not be expected that the data would lie on a continuous 
curve, unless the flow were either frozen or in equilibrium (without 
dissociation) . The bow wave shapes predicted from the theory and shown 
in figure 9 indicate that all of the bow waves for nonequilibrium flow 
become coincident with the wave for frozen flow as the apex is approached. 
Such coincidence of the maximum values of the wave angles with the values 
corresponding to frozen flow is not generally evident in the experimental 
data, primarily because accurate measurement of the wave angles could not 
be made from the shadowgraphs in the region near the apex (see fig. l^(b)). 

In most of the tests, the effect of dissociation was slight, and even 
at the highest temperatures (i.e., at the highest velocities), dissociation 
would produce less than a third of the possible variation of 0 V between 
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frozen and equilibrium at the available range of pressures of the tests: 
60 to 5*7 atmospheres behind the shock wave. It became evident that this 
limited amount of data applicable to oxygen dissociation would be inade- 
quate for the determination of the relaxation time for this process. 


Experimental Relaxation Time 


Reference 5 presents a discussion of the calculation of vibrational 
relaxation time from the theory of Landau and Teller (ref. 18 ) based on 
general quantum considerations. It is shown that Tp, the product of the 
relaxation time and the pressure of the relaxing gas, is a function of 
the temperature and of certain characteristics of the gases, including 
the collision cross sections and the characteristic temperature for 
vibration. For any one gas then, the theory shows that the relaxation 
time can be presented as a single curve, rp as a function of temperature. 
When the results of the tests are presented and compared with shock-tube 
results in this report, the product Tp is therefore used as the relax- 
ation time parameter. In figure l6 it is plotted as a function of the 
average temperature defined as follows, 

m _ T e + Tf 
2 

where T e is the calculated equilibrium temper at lire behind the shock 
wave and Tp is the temperature calculated assuming frozen flow in which 
the shock wave angle is the same as the angle in the shadowgraph from 
which the relaxation time is determined. The relaxation time parameters 
were obtained from the shadowgraph data as described in the section on 
reduction of the data. The circular symbols denote the tests in which 
the models were fired into still air and the triangular symbols denote 
tests with the models fired into the supersonic air stream of Mach 
number 3 - Symbols joined by vertical lines are from single shadowgraphs 
which yielded different relaxation times depending upon the radial position 
along the bow wave that was selected. The data shown were measured at 
three such positions, at the distances 0.6r^, 0.8r^, and r^ from the apex 

outward along the bow wave. Table I provides additional information (such 
as the local pressures and the equilibrium temperatures) pertaining to the 
data in figure l6. 

The triangular symbols in figure l6 (representing data obtained when 
the models were fired into the supersonic air stream) shown at temperatures 
of 1785° K, i860 0 K, 27 ^ 0 ° K, and 2830° K indicate somewhat longer relaxa- 
tion times than those predicted, corresponding to larger angles of the bow 
shock wave than those calculated from the shock tube relaxation time data. 

A similar result is indicated in table I and figure 15 (b) at a Mach number 
of 6.79 (and a temperature of 83O 0 K) for which the recorded bow shock 
wave angles were everywhere greater than the cone angle for frozen flow. 
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Apparently, when the models were fired into the supersonic air stream, 
the data are influenced hy factors other than just the relaxation time. 

Data shown in figure l 6 and table I at four temperatures (i860 0 , 

2220°, 3550°, and 4030° K) indicate that the value of rp varied 
considerably, depending on the choice of the portion of the bow wave 
upon which the calculations were based. This large variation corresponds 
to a large amount of curvature of the bow wave, the cause of which could 
not be ascertained but which may be partly a result of large yaw angles 
(up to about 4.4°) that were recorded when some of these shock waves were 
photographed, and partly a result of inaccuracy due to poor definition of 
the shock wave in the shadowgraphs. Contributing to the lack of definition 
were: occasional underexposed pictures resulting in poor photographic 

contrast; occasional fogging of the film plates by the gun muzzle flash; 
finite light duration causing a blurred shock wave image; finite size of 
the spark source which reduced image sharpness. 

The remainder of the data obtained with the models launched into 
still air and one set of data obtained using the supersonic air stream 
( shown in figure 17 at a temperature of 2090° K) , indicate reasonably 
consistent values of Tp as a function temperature . 


Comparison With Results of Other Tests 

Figure l6 also shows the vibrational relaxation times for pure oxygen 
and pure nitrogen measured in the shock tube tests described in refer- 
ence 5 • It is seen that for nitrogen rp is greater than for oxygen by 
more than an order of magnitude. It would be expected from this that the 
relaxation of air would not occur as a single process, but instead might 
be two processes: first, the excitation of the oxygen and an approach to 

a quasi-equilibrium state in which the nitrogen is almost frozen, and then 
as the excitation of the nitrogen progressed, an approach more slowly to 
complete equilibrium. If the processes actually took place separately in 
this way and if the flow distances were of the proper length, the bow wave 
shapes would show separate effects which could be identified with each of 
the two processes. Attempts to identify such effects in the shadowgraphs 
were not successful, however. The values of rp shown in figure l6 for 
the tests described in this report were calculated as if the relaxation 
were a single process, and, therefore, rp is an effective or apparent 
relaxation time parameter which characterizes the combined process in air. 
Since a comparison with shock tube results was considered one of the objec- 
tives of the experiments, similar "apparent" relaxation times were calcu- 
lated from the shock tube data. They are referred to in the following 
discussion as the predicted relaxation times. 

The predicted effective relaxation times in air have been calculated 
by two methods; in both the air is considered to be a mixture of oxygen 
(21.4 percent) and nitrogen ( 78 . 6 ), and the assumption is made that the 
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variation of the state of the gas mixture in the relaxing process is the 
sum of the variations of the states of component gases, as if they were 
not mixed. In the first method, which does not depend on the gas flow, 
the calculated relaxation time is defined as the time required for the 
mixture to vary its temperature (and density) from frozen to a specified 
fraction of the total change between frozen and equilibrium. In the 
second method, the bow shock wave is constructed and an apparent relaxation 
time is calculated from this shock wave. 


In using the first method, it is assumed that the density variations 
of the oxygen and the nitrogen can be expressed separately as exponential 
functions of time, so that the variable density of the mixture is given 
by the following equation. 


/ PM ~ P jj \ 

VV P V 







(16) 


where the subscript M refers to the gas mixture, 0 refers to oxygen 
and N to nitrogen. 'When t is equal to the relaxation time 

M 

(p - pf* ) is equal to l/e(p e - pf* ) , according to the definition of t 
M M MM 

that has been mentioned previously. Since the mass fractions of the two 
component gases are constant, so that 

P = 0.21k P P = 0.786 P 

0 e M e N e M 


p f = 0.214 p , p = 0.786 p 
0 


equation (l6) may be written 


0.214 ^e T M P//t 0 P ^ + 0.786 ^e T M P//t N p ^ = 0.368 


(17) 


If values for tqP and r^p are taken from the figure 16 (for oxygen 
vibration and nitrogen vibration, respectively) , noting that TqP is 
approximately equal to 0.04 r^p, equation (17) reduces to the relation 


TjvjP = 0.759 t n p ( 18) 

This is shown as the dashed curve in figure l6. Essentially the same 
relation as equation (l8) is obtained if the relaxation time of the 
oxygen is assumed to be zero. 
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The predicted relaxation time for air, calculated by the second 
method is described in detail in appendix D. The shape of a bow shock 
wave is calculated and from this wave shape, the relaxation time t and 
the parameter Tp are determined in the same way as in the case of the 
shock waves obtained experimentally. Results of the calculations are 
presented in table I and are shown as the rectangular symbols in figure 1 6 . 
They may be compared with the experimental results shown in this figure at 
the same temperatures. The variation indicated by the elongation of these 
rectangular symbols results from calculating the effective rp for several 
points along the wave. As the location of the point on the shock wave 
varies, the relative contribution of the oxygen and the nitrogen relaxation 
processes changes so as to cause a variation in the effective relaxation 
time parameter. 

The shock-tube tests reported in reference 5 provide data for 
vibrational relaxation time in oxygen only up to a temperature T of 
3000° K. The results were extended to the higher temperatures indicated 
in figure 1 6 by plotting the log of Tp as a function of T -1 / 3 (as 
suggested in ref. 5 ) and extrapolating. At these high temperatures, the 
predicted values of rp are quite insensitive to variations of the 
relaxation time of the oxygen. 

The apparent values of rp obtained from the tests would be 
expected to have the same characteristic as the values predicted by the 
method involving the construction of a bow shock wave, that is, a variation 
with the location of the points selected on the shock wave. In addition, 
these values are also to some extent functions of the pressure and the 
magnitude of the characteristic flow dimension (i.e., the model size). 

It appears that for the majority of the points shown, the values of the 
relaxation time parameter determined experimentally agree fairly well with 
values calculated by either of the two methods described, in which the 
results of shock-tube tests in pure oxygen and pure nitrogen are combined 
to furnish values for relaxation times in air. 

In common with most experimental techniques that have been employed 
to measure relaxation times at high enthalpies, the results obtained in 
the tests described in this report do not provide precise numerical values 
for Tp* However, in view of the very large differences in the values of 
this parameter that have resulted from use of various theories for calcu- 
lating relaxation times for some of the processes involving adjustment 
rates of the various internal degrees of freedom of a gas, in particular, 
for dissociation and the other chemical reactions, an experimental method 
which yields values of Tp even within an order of magnitude can be of 
considerable value. The results shown in figure 1 6 indicate a mean varia- 
tion of Tp with temperature from which numerical values oan be specified 
well within an order of magnitude. 

The tests reported in reference 5 indicated that the vibrational 
relaxation time of oxygen in a nitrogen mixture was increased because the 
oxygen-nitrogen collisions are only about 4o percent as effective in 
transferring energy as are the oxygen-oxygen collisions. If this factor 
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is included in calculating the effective relaxation times of air, the 
values of t for oxygen in air would he about twice those shown in 
figure 1 6 for pure oxygen. All of the effective values of t calculated 
from the shock-tube results would be increased slightly, but the general 
comparison of such values with the values obtained from the tests reported 
herein would not be significantly different from that shown in figure 1 6 , 
because oxygen is a minor constituent (21 percent) of air. 

Results of the tests in still air which covered a range of tempera- 
tures (T) from 1990° to ^260° K and velocities from 7^260 to 12,000 feet 
per second, indicated values of rp from 60.7 to O .78 microsecond- 
atmospheres, corresponding to a range of relaxation times from b.^b to 
0.16 microseconds. At the highest velocities, about 30 percent of the 
change in bow wave angle between frozen flow and equilibrium is caused by 
oxygen dissociation, and the relaxation time for the dissociation would 
have some influence on the effective relaxation time for air calculated 
from the experimental data. The results thus far obtained, however, are 
not sufficiently precise to be used to calculate a relaxation time for 
oxygen dissociation. 


Effect of Water Vapor 


The effective relaxation time for air is undoubtedly influenced by 
the presence of gases other than oxygen and nitrogen. Reference k indi- 
cates that water vapor in ambient air can have an important effect on the 
relaxation time for oxygen vibrations, because of the much greater effec- 
tiveness of oxygen-water molecular collisions compared with oxygen-oxygen 
and oxygen- nitrogen collisions . The absolute humidity of the air in which 
models were flown was always about equal to that of ambient air, when the 
tests were conducted in still air. In these tests it was not feasible to 
use dry air because there was always some leakage of ambient air into the 
test section whenever the test section was partially evacuated. It is 
indicated in reference 5 that an amount of water vapor up to 3 parts in a 
thousand did not have a measurable effect on the relaxation time in oxygen. 
In the present tests the amount of water vapor was as much as four times 
that in the tests reported in reference 5 and may have had some effect. 

The data have been examined for any consistent relation between the 
absolute humidity of the ambient air and differences between the predicted 
and experimentally determined relaxation times. No consistent correlation 
was evident. 


CONCLUSIONS 


Conical models of low-fineness- ratio have been tested in free flight 
at velocities from 4,600 feet per second to 12,000 feet per second, in 
which the calculated equilibrium temperatures behind the shock wave ranged 
from 1510° R (8Uo° K) to 6200° R (3^0° K) and temperatures in frozen flow 
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ranged from 1580° R (880° K) to 9l4o° R (5070° K) . Analysis of the 
results of these tests has led to the following conclusions: 

1. The shape of the how shock wave can he used to calculate a 
relaxation time in a nonequilibrium flow in the heated region behind the 
shock wave. 

2. Comparison of the relaxation times in air calculated from the 
present tests with effective relaxation times estimated from results of 
shock-tube tests with pure oxygen and pure nitrogen indicated that in the 
temperature range where vibrational relaxation would be expected, these 
two experimental techniques yield results that agree within the experi- 
mental scatter. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif May 9, i960 
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APPENDIX A 

f 

DEVELOPMENT OF APPROXIMATE THEORY FOR CALCULATING 
EQUILIBRIUM FLOW ON BLUNT CONES 


The basic equations of axisymmetric cone flow from reference 11 may- 
be written, with the notation indicated in figure 2(a), 


du 

d0 


v = 0 


(Al) 


d v 
d 0 


+ u 


+ 


1_ dp _ 
pv d0 


(A2) 


(pv sin 0) + 2pu sin 0=0 

U.C7 


(A3) 


Since flow along streamlines is isentropic behind the shock wave, the 
sonic velocity given by the relation a 2 = (dp/dp)^, may be expressed in 
the relation 


dp _ 2 dp 

de de 


(Ah) 


The following expression is obtained by combining the above relations. 


du 

d0 


dfu 

d0 2 


+ u = 


dfu 

d0 2 


+ | cot 0 + 


2u 


(A5) 


At the solid cone surface, designated by the subscript s, v s = 0, and 
so from equation (Al), (du/d0) g = 0. Equation (A5) at the solid cone 

surface then becomes 


($).- (A6) 

Expanding v in a series referred to its value on the cone surface, and 
dropping the higher order terms (assuming d 2 u/d0 2 constant) 



(e - e s ) 


V 


s 


(A7) 


2k 




From equations (A 6) and (A7) 

v = -2u s (e - 0 S ) (A8) 

Comparison of values obtained from this approximate relation with 
values given by reference 11 indicates good agreement for the blunt cones 
considered here. It was noted that even better agreement results if the 
local radial velocity u is used in equation (A8) instead of the velocity 
Ug on the surface; that is, defining A0 = 0 V - 0 S 


v w = -2u„A0 (A9) 

If the velocity ahead of the shock wave, U x , is resolved into its 
normal component u ni and tangential component u t x > an< ^ ‘the velocity 
behind the shock wave into components u^ and u-^, equation (A9), 
together with oblique wave relations, provides the following: 

Since = u^ and v w = -u^ = -2u„A0 

u n2 = 2u t2 A0 (A10) 


tan 0 W 


Un i u n x 


(All) 


If a 

wave. 


is the angle of the streamline deflection through the shock 


tan (0 W - a) 


^2 


(A12) 


2u t2 A0 

U t 2 


2 ( 6 w " 6 s> 


(A13) 


from equations (A10) and (A12) . Also, from equations (All) and (A12) 


Ur- 


tan (0 W - a) = tan 0 W 


•n. 


Then, from equations (A13) and (AlU) 


(Al4) 



tan 0 W = 2(0 W - 0 S ) 


(A15) 


-o ro ro > 


l\J l\J 


4Z 


25 


Inserting the relation from continuity requirements 


Un 2 Pi 
u n, P, 


equation (A15) becomes 


py ) "* :an ®w - 2 (0 W — 0 S ) 


(Al6) 


Equation (Al6) relates two of the variables we are most interested in, 
the wave angle 0 W and. the density ratio p 1 /p 2 - Since both are unknown 

we require another relation between them and for this purpose write the 
momentum and energy equations for flow through the oblique conical shock 
wave using the normal component of the Mach number M n and assuming the 
gas in the free stream is ideal, 


rT = 1 + ( 1 - 7^ ) /Mn' 


sr 1 + 


1 - 


7M n 2 

(A17) 

) 

(A18) 


The equation of state may be written 


tp 


(A19) 


where t = Tg/l^ and Z 2 is the compressibility, P 2 /p 2 RT 2 ' Equations 
(A17), (Al8), and (A19) may be combined to give the relation 


+ £1 

^7 ( Rt l ) ( f ^ 

\P2 / P 2 

y I ± ) \ XZj 2 


- tz 2 = 0 


(A20) 


where (3 = (h 2 /h 1 ) (T^^) . If an effective ratio of specific heats, 7 E > 
is defined by the relation 


^E + 1 
?E- 1 


7 + 1 t 
7-1 + t - 1 


(P - 1) 
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7-1 


- (Z 2 - 1) 1 + 


(A21) 
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equation (A 20 ) becomes 


+ r 7 E + 1 

?*) p 2 l^e - i 


(t - 1) 


- t = 0 


(A 22 ) 


It may be noted that in the special case of a perfect gas (p = 1, Z 2 = l) , 
the definition of 7g, equation (A 2 l), reduces to 7g = 7. Neither P 

nor Z appears in equation (A 22 ) , which is the same as the expression 
for the density ratio in a perfect gas having 7 = 7^, and the same value 

Jli 

of t . 


An expression for the pressure is found by writing the momentum 
equation for the normal component across the shock wave 


P, - P~ = P,u n p 


2\ 


From continuity. 


PiX 


= P p u n„ 


so that 

P 1 ~ P 2 _ _ Pi 

P l Un l 2 " P 2 

Since, by definition, 

_ P 2 -Px 
P (l/ 2 )p 1 U 1 2 


and 


un x = U x sin 


The pressure coefficient may be expressed as follows: 


C 


P 


2 sin 2 0 w 



( p i/ p 2 ) 


(A 23 ) 


In the hypothetical flow, the enthalpy ratio is equal to the 
temperature ratio. Putting t = h 2 /h. and M n = M, sin 0 W into 

equation (Al8), 
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M lE 2 sin 2 0 w = 


t - 1 


1 - (p x /p 2 ) £ 


(1/2) (r E - i) 


Writing equation (Al 8 ) for the general gas with h z /ln. x = 0t, 


M 1 2 sin 2 e w = 


Pt - 1 


1 ~ (P x /P 2 ^ 


( 1 / 2) ( 7 - 1 ) 


The Mach numbers of the two flows are related by the following 
obtained by combining equations (A 2 ^) and (A25) • 


Mi 



" (Pt - 1) (y E - 1) " 
_ (t - l) (7 - l)_ 


2 


(A24) 


(A25) 


equation. 


(A26) 
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APPENDIX B 
EXAMPLE CALCULATION 


This appendix presents an example of the calculation of the bow 
shock wave angle and the pressure ratio across the shock for a 

52. 5° cone. The temperature and pressure behind the shock wave are 
assumed to be, respectively, 38OO 0 K (6840° R) and 10 atmospheres, and a 
free-stream temperature of 300° K (5^0° R) is selected. The enthalpies 
and compressibility factor for air obtained from references 1 and 3 were 
as follows for the assumed conditions. 

h 1 = 3.83 RT 0 

h 2 = 74.1 RT 0 
Z 2 = 1.066 

From equation (A20) , p /p 2 = 0.1179. From equation (A2l) , assuming 
7 = 1 . 405 , the value 7^ = 1.244 is obtained. For this value of 
7 e and 0.3t = 0.3(T 2 / t i) = 0.3(3800/300) = 3.8, figure 3 gives the 
the effective Mach number 

m e = 11.62 

With this M E , equation (A26) gives the stream Mach number 

Mi = 11.29 

From figure 4 and the above values of 7^ and 

A0 = 5*50 


The bow shock wave angle is 


6 w = 0 S + A0 
= 52-5 + 5.50 
= 58.00 

Obtaining the pressure coefficient from figure 7, 


Cp = 1.266 
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This together with the assumed pressure behind the shock wave permits 
the calculation of a pressure ratio and a stream static pressure 

CpM^r 

2 

(1.266) (11. 29) 2 (l.409) 

2 

= 113. 4 

p i = P 2 AP 2 /P 1 ) = 10/113.4 

= 0 . 0882 atmosphere 
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APPENDIX C 


EFFECT OF THE BOUNDARY LAYER ON THE SHOCK-WAVE SHAPE 


The magnitude of the change in the shock-wave shape due to the 
boundary layer would depend on the condition of the boundary layer. In 
some of the tests, it was possible to determine whether the boundary 
layer was laminar or turbulent by inspection of the shadowgraphs . Such 
tests indicated that in most cases the boundary layer was entirely laminar 
however, turbulent boundary layers were observed in a few cases, and the 
possibility existed that the boundary layers were turbulent in some of the 
tests for which it was not possible to determine their condition from the 
shadowgraphs. In order to establish the limits of the possible effects, 
the maximum effects of laminar and turbulent boundary layers on the shock 
wave shapes were calculated. 

From heat-transfer calculations (and also from experimental results) 
it can be shown that the model surface temperatures in tests such as those 
discussed in this report are always considerably lower than the adiabatic 
wall temperature, because of the short duration of the tests and the high 
conductivity of the model. In calculating the thickness of the laminar 
boundary layer, the effect of heat transfer to the model surface has been 
taken into account by use of the results presented in reference 19. In 
this reference, theoretical results are given for an assumed ratio of wall 
temperature to temperature outside the boundary layer of 0.25 . In most of 
the tests, the temperature ratio was estimated to be even smaller than 
this, but this ratio is typical of the test conditions for which the 
boundary-layer thickness was maximum. In addition to assuming that the 
temperature ratio of 0.25 would be applicable, a combination of conditions 
was selected in which other factors would cause the boundary-layer thick- 
ness to be large, and the boundary-layer thickness was calculated for this 
one case to determine the maximum effect of the laminar boundary layer. 

For this combination of conditions, the Reynolds number behind the bow 
shock wave was 0.9 million per inch, or a Reynolds number of one million, 
based on the slant length of the cone. (The Reynolds number based on the 
stream conditions was 1.7 million per inch.) 

The thickness of the boundary layer on a cone was related to that on 
a flat plate by the relation 


which results from the rule (given, e.g., in ref. 20) that the skin 


boundary-layer displacement thickness at the most rearward location on 



friction on a cone is 


times the value for a flat plate. The 
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the conical surface was calculated to be 0.00053 inch, corresponding to 
an effective angular change in the model cone half angle of 0.025 . This 
change in the effective model cross section would produce a negligible 
change in the shape of the bow shock wave. Although the way in which the 
effect of heat transfer is taken into account is not exact insofar as the 
surface temperature of the model is not known, if the effect were entirely 
disregarded, the calculated thickness of the laminar boundary layer would 
be no more than twice the value obtained when wall cooling was assumed 
and might still be neglected within the accuracy of the measurements. 

To estimate the maximum effect of a turbulent boundary layer, as in 
the case of the laminar boundary layer, calculations were made for one 
combination of conditions in which the displacement thickness might be 
expected to be a maximum. The same Reynolds number was selected as in 
the case of the laminar boundary layer. The thickness of a turbulent 
boundary layer on a cone was calculated by an unpublished method in which 
the effect of heat transfer can be taken into account if the effect of 
heat transfer on the skin- friction coefficient is known. The skin friction 
for flow with heat transfer was computed by the method described in refer- 
ence 21. For the conditions considered, the maximum calculated displace- 
ment thickness was 0.0016 inch. The effective increase in the cross 
section of the model cone corresponds to an increase in the solid cone 
angle (in the region where shadowgraph data were measured) of about 0.12° . 
The influence of this incremental cone angle on the bow shock wave angle 
is indicated in figure 17, which shows the bow wave angle as a function 
of cone angle for frozen flow at Mach numbers of 6, 10, and infinity, and 
for equilibrium flow at Mach numbers of 6 and 10. The maximum slope of 
these curves, for the range of conditions of the tests, was 2.2. This 
indicates that the angular coordinate of the bow wave could increase by 
0.26 , if the model has a turbulent boundary layer originating at the apex 
of the cone. This increase would have some effect, such as to increase 
the calculated value of the relaxation time, that would increase with a 
decreasing difference between A 0^ for frozen flow and A 0 e for 
equilibrium flow. The largest effect then occurs at the lower model 
velocities, and for the smaller solid cone angles. 

.o 

Because this value of the angular increment, 0.26 , results from a 
combination of assumed conditions such as to produce the maximum effect, 
and because it is believed that in most of the tests the boundary layer 
was partly or entirely laminar so that the effect would be much smaller, 
no correction was made for the effect of the boundary layer on the shape 
of the bow shock wave. 
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APPENDIX D 

PREDICTED EFFECTIVE VIBRATIONAL RELAXATION TIMES 
BASED ON A CALCULATED SHOCK-WAVE SHAPE 


The second of the two methods employed in this report to calculate 
the predicted relaxation times for air is presented in this appendix. 
First, relations are obtained from which the bow wave shape can be 
calculated, and then the relaxation times are calculated from this wave 
shape. Calculations are made for conditions selected so that the results 
can be compared directly with each set of experimental results. 


The following equation is assumed to specify the time variation of 
temperature of the air as it approaches equilibrium after having been 
heated as an inert gas by the shock wave. 


T - T e = (Tf - T g )e 


-t/' 


0 


(T g - T e )e 


-t/r, 


N 


(Dl) 


where the subscript g refers to the state of the gas corresponding to 
oxygen in equilibrium while the nitrogen remains frozen. Introducing the 
density and pressure from the equation of state into equation (Dl) 


JL - Is. = (SL - e' t/T ° e‘ t/TB 

pR p e R \Pf R PgR/ \Pg R Pe R / 


(D2) 


At any given free-stream Mach number, it can be shown that the pressure 
ahead of a blunt body is not a sensitive function of the gas state, in 
which case it can be assumed that p, p e , and p^ are equal. Also, with 
no dissociation, the gas constants are equal. Equation (D2) can then be 
written 





(D3) 


It is now desirable to relate the bow shock wave geometry to the 
density. The shock wave angle is given as a function of density in 
equilibrium flow by equation (Al6) , which can be written 


_1 = _2_ A0 
P 2 P -l tan 9 W 


W 


ro ro > 
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Putting equation (D4) into equation (D3) and observing that p x , the free 
stream density, is the same in each term, 

AO A0 a f A0f A0g \ "^a/ r Q / A0g AOq \ ”^a/ T ]^ 

tan 6y tan 6 e \tan Of tan OgJ \tan 0g tan 6 e J 6 

(D5) 

For this expression to have a meaning, it is necessary to assume that an 
average effective density can be defined which is related by equation (d 4) 
to the angular coordinates of points on the shock wave in nonequilibrium 
flow. This density is, in general, a function of the radial coordinate. 

In equation (D5) the time t a is the time required for a small volume of 
gas to attain the average effective density as it moves downstream from 
the shock wave. 

To determine whether it is possible to assume that an effective 
average density and a corresponding flow time t a can actually be 
defined, equation (DU) was applied to a nonequilibrium flow, assuming 
that the gas undergoes only a single relaxation toward equilibrium. The 
result was compared with the result obtained using the derived relations 
for nonequilibrium flow and figure 9* The purpose of this comparison 
was to see if any value of t a would bring the results into agreement. 

It was found that relatively good agreement is obtained (at least for the 
case of model cone angles of 52 . 5 ° and 55°) if t a is taken as 0.2(i^/uq). 

Equation (D5) 'was then used with this value of t a to calculate 
A0/tan Qy for three radial locations on the shock wave, 0.6^, 0.8^, 
and r^. For each of these values of A0/tan 0 W , the angular coordinate 
Gy was calculated and y/y^ was computed: 

y _ ^w ~ 0 e r 
~ 0 f “ 0 e r b 

The relaxation time ratio r^/ru Q was determined from figure 9 and 
converted into the parameter Tp using the same values of r ^ , u Q , and p 
as those used in reducing the corresponding experimental results. 
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TABLE I.- TEST CONDITIONS AND RELAXATION TIME PARAMETERS 



Mi 

* x > 

°R 

V 

°R 

Tj 

°R 

T, 

°K 

Pressure 
behind shock 
P 2 , atm 

Average 
velocity 
behind shock 
Uq, ft/sec 

e B’ 

deg 

Model 

reference 

length 

ft 

Radial 

coordinate 

r / r b 

TP* 

experiment 

^sec-atm 

TP* 

predicted 
u sec- atm 

Tp( experimental) 

Tp(predicted) 

Still 

6.48 

522 

3360 

3590 
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38.9 

3610 

52.7 

0.0762 

0.6 

48.3 

33.6 

1.43 

air 










.8 

41.7 

27.9 

1.49 
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1.0 

42.0 

14.0 
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3500 
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O Exact solution, # s =50° y - 1.405 
□ Exact solution, # s =55°, y- 1.400 

Approximate solution, 9 S Z 50°, y- 1.405 

Approximate solution, # s = 55°, y- 1.400 
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Figure 1.- Comparison of the calculated bow shock wave angles from approximate theory with the 

wave angles from exact theory; frozen flow. 
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(b) Nonequilibrium flow. 

Figure 2.- Sketch showing model geometry and flow field notation. 
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Figure 3 . 
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Graph showing relation between temperature ratio t, effective y , 
and effective Mach number; 0 S = 52.5° • 
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Figure 4.- Graph relating incremental shock wave angle, effective 7 , and 

effective Mach number; 0 S = 52.5° • 
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Figure 5 



- Graph relating temperature ratio t, effective 7, and effective 
Mach number; 0 S = 55° • 
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Figure 6 .- Graph relating incremental shock wave angle, effective 7 , and 

effective Mach number; 0 g = 55° • 
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Figure 7.- Pressure coefficient behind the bow shock wave. 
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Figure 8.- Calculated incremental bow shock wave coordinate angles for frozen and equilibrium flow 

in air; 9 S = 52.5°, T x = 5^0° R. 
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Figure 9*" The distance between the bow shock wave in nonequilibrium flow 

and the equilibrium- flow shock wave. 
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Figure 10.- Photograph of three cone-cylinder models used in free- flight tests. 
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Figure 11.- The ratio of the yaw angle of the conical shock axis to the yaw 

angle of the solid cone; y = 1.405* 
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Figure 12.- Increase in apparent cone angle of the bow shock wave as a result of yaw angle of the 

conical shock wave axis . 
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M = 6.53- 


M = 7.24. 


(a) Still air tests. 
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(b) Supersonic air stream; M = 10-99* 

Figure 13*- Typical shadowgraphs of cone-cylinder models in free flight 
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(a) Rectangular coordinates of shock wave. 
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(b) Polar coordinates of shock wave. 

Figure l4.- Example of bow shock wave plotted from shadowgraph data; 

0 S = 52. 5°, M x = 6.53. 
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(a) e s = 52. 5°; T x = 5^0° R. 


Figure 15-- Comparison of the experimental and theoretical incremental how shock wave coordinate 

angles for frozen and equilibrium flow. 
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(b) e s = 52 . 5 0 ; t 1 = 193° R 

Figure 15.- Continued. 
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Shock wave angle relative to Shock wave angle relative to 

cone surface, A deg cone surface, A#, deg 
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(c) e s = 55°; t x = 5^0° r 
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Frozen flow, y - 1 .405 

Only oxygen vibrations excited, no dissociation 

Oxygen and nitrogen vibrations excited, no 

dissociation 

Equilibrium 

O Experimental results 
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(d) e s = 55°; i 1 = 193° R 


Figure 15.- Concluded. 
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Figure 16.- Apparent relaxation time parameter as a function of average 

temperature behind the shock wave . 
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Figure IT Variation of bow shock wave angle with solid cone angle for 
various Mach numbers and free-stream temperatures. 


VJ1 

VJ1 


